Nuclear shadowing and extraction of Ff — F£ at small x from deuteron collider data 
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We demonstrate that leading twist nuclear shadowing leads to large corrections for the extraction 
of the neutron structure function F£ from the future deuteron collider data both in the inclusive 
and in the tagged structure function modes. We suggest several strategies to address the extraction 
of F£ and to measure at the same time the effect of nuclear shadowing via the measurement of the 
distortion of the proton spectator spectrum in the semi-inclusive eD — > e'pX process. 
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One of the important components of the electron-ion 
collider project (EIC) 

is the measurement of the non- 
singlet structure function of the nucleon, F% — down 
to rather small values of Bjorken x. A similar program 
is discussed for HERA III. The main objective of us- 
ing deuteron beams at HERA is obtaining more accurate 
information on the nonsinglet structure function of the 
nucleon at small x. Since Fl^ and F£ are rather close at 
small x, nuclear shadowing significantly affects F% — i 7 ^ 1 , 
when extracted from the deuterium data. Indeed, while 
R = (Ff - F?)/F% < 4 x 1CT 2 for x < 1(T 2 , the shad- 
owing correction modifies by 1-2% and, hence, R is 
modified by twice as much. Therefore it is crucial to de- 
termine the value of the nuclear shadowing correction, as 
well as its uncertainties, both to the inclusive structure 
function F® and to the tagged structure function, when 
the spectator proton (a proton with momentum < 0.1 
GeV/c in the deuteron rest frame) is detected ensuring 
the kinematics maximally close to the scattering off a free 
nucleon. 

In this paper, we analyse both of the above mentioned 
structure functions using the theory of leading twist nu- 
clear shadowing Q , which is based on the existence of the 
deep connection between the phenomena of high energy 
diffraction and nuclear shadowing demonstrated within 
the reggeon calculus for hadron-deuteron scattering by 
Gribov |3] , and the Collins factorization theorem for hard 
diffraction . In the case of the tagged structure func- 
tions, we also use the Abramovskii, Gribov, Kanchcli 
(AGK) cutting rules 5], which relate shadowing effects 
for the total and partial cross sections. We formulate the 
requirements necessary for the extraction of F% from the 
deuterium data and demonstrate that a detector with 
good proton momentum resolution will be able to deter- 
mine the value of the shadowing correction. 

We begin by discussing nuclear shadowing in inclusive 
DIS off deuterium. On the qualitative level, the Gribov 
result for nuclear shadowing off the deuteron |3| could be 



understood as a consequence of the interference between 
the amplitudes for diffractive scattering of the projectile 
off the proton and off the neutron of the deuterium target. 
Such interference is possible for small x, x < 5 x 10~ 2 , 
where the minimum momentum transfer to the nucleon, 
~ xttin, becomes smaller than the average nucleon mo- 
mentum in the deuteron. The corresponding double scat- 
tering diagram for the "f*D scattering is presented in 
Fig. ^ The nuclear shadowing correction to the total 
j*D cross section is given by the imaginary part of this 
diagram, which is obtained by making all possible unitary 
cuts (denoted by dashed lines in Fig. [I}. The application 
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FIG. 1: All possible unitary cuts of the interference diagram 
giving rise to nuclear shadowing correction to 7*D cross sec- 
tion. 



of the AGK cutting rules explicitly demonstrates that 
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the imaginary part of the interference graph decreases 
the total 7*£> cross section by the factor proportional to 
(1 — ?7 2 )(Im.4) 2 , where A is the amplitude for the photon- 
nucleon diffractive scattering and r\ = KeA/lmA. 

A consideration of the individual unitary cuts in Fig. 2] 
demonstrates that the interference diagram decreases the 
cross section of inelastic interactions with a single nu- 
cleon (top right and bottom left graphs) by the factor 
proportional to 4(Im.4) 2 , and also results in simultane- 
ous inelastic interactions with two nucleons (botton right 
graph) with UdouUe = 2(Im_4) 2 . Also, the interference di- 
agram increases the probability of diffraction on the nu- 
cleus by the factor proportional to (1 + ?7 2 )(Im.4) 2 (top 
left graph), see Ref. 0. 

Within the framework of the Gribov theory of nuclear 
shadowing, the deuterium inclusive structure function 
reads 



F 2 D (x,Q 2 ) = FP(x,Q 2 )+FZ(x,Q 2 ) 



xpoiMt +4(x JP TO Ar ) 2 ) 



dx P dq 2 F° {A) ((3,Q 2 ,x P ,t) 



(1) 



where F D ^ is the nucleon diffractive structure function, 
Pd is the deuteron form factor, and \t\ = q 2 + (xpuin) 2 ■ 
Since the t-dependence of po is rather moderate (com- 
pared to heavier nuclei) , the integral in Eq. Q is sensitive 
to F D{i \t) up to -t < 0.05 GeV 2 . Also, in Eq. Q, the 
factor (1 — ?7 2 )/(l + ?7 2 ), where 
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2 d\n{l/x P ) 



: (a P (t = 0)-l) 



(2) 



accounts for the real part of the amplitude for the diffrac- 
tive scattering and leads to a reduction of the shad- 
owing by about 20%. Note that in the reggeon calculus 
derivation it was assumed that r\ = 0, which is natural 
for the amplitudes slowly increasing with energy. This is 
not the case for DIS and, hence, r\ should be taken into 
account. 

Using the QCD factorization theorem for hard diffrac- 
tion [j], it is possible to extend the Gribov theory in 
order to calculate nuclear shadowing for the quark and 
gluon parton densities of the deuteron, fj/D, at small x 
as follows, see Ref. Q and subsequent publications 0, , 

fj/D{x,Q 2 ) = fj/ p (x,Q 2 ) + fj/ n (x,Q 2 ) 
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/ dxpdq 2 fjj N (f3,Q 2 ,xp, t) 

J X 



xp D (4q 2 +A{xpm N ) 2 ) 



(3) 

The results of the calculation of the ratios Ff / (F% + F£) 
and g D /(2g N ) using the HI diffractive fit for F D{ ?> 
are presented in Figs. [21 and [3] for a range of x and Q. In 
this calculation we used the Paris NN potential. One can 




FIG. 2: The ratios F 2 D /(F 2 P + F 2 ") and g D /(2g N ) as functions 
of x. The solid curve corresponds to Q = 2 GeV; the dashed 
curve corresponds to Q — 5 GeV; the dash-dotted curve cor- 
responds to Q = 10 GeV. 



see that substantial shadowing is expected for the small 
x region. 

Note also that the nuclear shadowing correction to F^ 
for Q 2 < 1 GeV 2 is expected to be substantially larger 
than that for Q 2 ~ 4 GeV 2 . The enhancement of diffrac- 
tion at small Q 2 by higher twist effects such as, for in- 
stance, vector meson production, will increase nuclear 
shadowing up to a factor of two. Note that the applica- 
tion of the Gribov formalism to nuclear shadowing in the 
NMC kinematics (Q 2 < 2 GeV 2 ), where the double scat- 
tering term dominates, leads to a very good description 
of the data, see e.g. [11|. 

Since the diffractive cross sections are likely to be 
known with accuracy of about 10% for small t, it ap- 
pears that the accuracy of the calculation of the nuclear 
shadowing correction to the inclusive cross section will 
be better than 0.2% for a wide range of Q 2 . Correspond- 
ingly, the theoretical uncertainty for the ratio of F£ /Ff 
will not exceed 0.4%, which is likely to be smaller than 
possible experimental systematic errors. 

A strategy, which is complimentary to the inclusive 
measurement of F% , is the use the neutron and proton 
tagging. We will primarily focus on tagging of scattering 
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FIG. 3: The ratios /(F% + F?) and g D /{2g N ) as func- 
tions of Q. The solid curve corresponds to x — 10~°; the 
dashed curve corresponds to x — 10 -4 ; the dash-dotted curve 
corresponds to x = 1CP 3 . 



off a neutron via the detection of the spectator protons 
and will only briefly comment on other possibilities. 

It was pointed out above that according to the AGK 
cutting rules the overall correction of the nuclear 
shadowing effect for nondiffrative events is four times 
larger than for the inclusive scattering. However, it is 
concentrated at relatively large transverse momenta (fc t ) 
of the protons and depends on k t rather strongly. Hence, 
two strategies will be possible. One would be to select 
only very low pt protons with a gross loss of statistics. 
The other, more promising approach is to measure the 
Pt dependence of the spectrum up to pt ~ 250 MeV/c 
and then to use this for the verification of the theory 
and for the measurement of the effects of the rescatter- 
ing (nuclear shadowing) in an independent way with the 
subsequent correction of the data for this effect. Provided 
good momentum resolution of the proton spectrometer, 
one would be able to make longitudinal momentum cuts 
to suppress/increase the shadowing effect (the shadowing 
effects are minimal for \2p^ jpu — 1| < 0.1). 

The following analysis demonstrates the dependence 
of the nuclear shadowing correction to the "f*D — > pX 
cross section on the transverse momentum of the specta- 



tor proton p t . The impulse approximation expression for 
the differential cross section of interest reads 



(Pp 



IA 



l (2-2x L )(u 2 (p)+w 2 (p)), (4) 



where ij, = E p /Ejj w (l—p z /mN)/2 is the Feynman x of 
the spectator proton; p = (pt,Pz) is the three-momentum 
of the detected (spectator) proton in the deuteron rest 
frame; u and w are the S- wave and D-wave components 
of the deuteron momentum space wave function normal- 
ized as 47T J dpp 2 (u(p) 2 + w(p) 2 ) = 1; the factor 2 — 2xl 
is the Miiller flux factor. The presence of leading twist 
nuclear shadowing adds a nuclear shadowing correction 
to Eq. @, so we find for the complete cross section 

^ = <P "(2 - 2x L )(u 2 (p) + w 2 (p)) ±* 

u{p)u(p') 



(pp 

x 



dx 



TP 



((3,Q 2 ,xjp,t) 



-w(p)w(p )(- - -) 

2 p^p' z 2 



(5) 



where p' = p + q t + (xpmN)e z . The additional factor 
of 3 — rj 2 in front of the shadowing correction is a re- 
flection of the AGK cutting rules. In the considered 
case, the shadowing correction is given by the sum of 
the two top graphs in Fig.Q] The left graph corresponds 
to the enhancement of the diffractive processes by the 
factor of (1 — 7] 2 )(Im^l) 2 , while the right graph corre- 
sponds to screening from the single inelastic subprocesses 
and is proportional to — 4(Im.4) 2 . In Eq. ©, the factor 
u(p)u{p') + w{ P )w(p')(3(p-p') 2 /(2p 2 p' 2 ) - 1/2) is simply 
the unpolarized deuteron density matrix. Also, we ne- 
glected production of deuterons in the diffractive channel, 
which would somewhat increase the shadowing effect. 

As one can see from Eq. l|B*jl. nuclear shadowing sup- 
presses the spectrum of the produced protons. This effect 
can be quantified by considering the ratio R of the com- 
plete expression given by Eq. JSJ to the impulse approxi- 
mation expression given by Eq. (JIJ. Figure 0] depicts the 
ratio R as a function of Bjorken x for p t — (0, 100, 200) 
MeV/c and p z — (top panel) and for p t = and 
\p z \ = 100 MeV/c (bottom panel). The calculation is 
made at Q — 2 GeV. 

Two features of Fig. 0] are of interest and importance. 
Firstly, nuclear shadowing works to decrease the ratio 
R as pt increases. This is expected from the general 
picture of nuclear shadowing since large pt corresponds 
to small transverse distances between the two nucleons, 
which leads to their shadowing, as was first pointed out 
by Glauber ^3] an d later generalized by Gribov [3} • Sec- 
ondly, the suppression of R at large p t is strikingly large. 
This is a common feature of semi-exclusive reactions with 
nuclei. Indeed, at large pt, while the impulse approxima- 
tion term is suppressed by the nuclear wave function, the 
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FIG. 4: The suppression of the proton spectrum by the nu- 
clear shadowing correction. Top panel: The solid curve cor- 
responds to p t — 0; the dashed curve corresponds to pt = 100 
MeV; the dash-dotted curve corresponds to pt = 200 MeV. 
Bottom panel: p t = and p z = (0, -100, 4-100) MeV/c. 



rescattering term survives and gives the dominant con- 
tribution. An example of this effect in electron-deuteron 
interactions in the TJNAF kinematics can be found in 
Ref. H3. 

Note also that there is a non-spectator contribution to 
the nucleon spectrum, which originates predominantly 
from diffractive scattering off the proton and dominates 
at large pt, Pt > 300 MeV/c. As a result of the AGK can- 
cellations, the non-spectator contribution is given by the 
impulse approximation. The contribution has a broad 
Pt distribution, oc e~ Bpt with B ~ 7 GeV -2 , and can 
be subtracted using the measurements at, for instance, 
Pt > 400 MeV/c. It appears that it will be possible to use 
most of the spectator protons from most of the deuteron 
wave function in the tagged method. Hence the spectator 
tagging will not lead to a loss of statistics. The method 
will allow to introduce corrections for nuclear shadowing 
for small p t with high precision, which would lead to the- 
oretical errors in the determination of F^ 1 at the level of 
a fraction of percent. 

One can also use simultaneous tagging of protons 



and neutrons, when both neutron and proton are de- 
tected in the reactions j*D — > nX and — > pX. 
In this case, nuclear shadowing will cancel in the ratio 
a 1 D ^ nX ja 1 d ~*p x and the main errors in the measure- 
ment of F£ will be due to the determination of relative 
efficiencies of the proton and neutron taggers. 

One could also try to obtain the ratio Fg /F% from the 
comparison of the rate of the tagged proton scattering 
events with the neutron spectator to inclusive e D scat- 
tering. Such a strategy could also have certain merits as 
it avoids the issue of luminosity and does not require a 
leading proton spectrometer. The disadvantage of this 
strategy is the sensitivity to the nuclear shadowing ef- 
fects and errors in the acceptance of the neutron detec- 
tor. One possible way to deal with the latter problem will 
be to perform measurements at very small x and large 
energies, where the ep and en cross sections are equal 
to better than a fraction of 1% and, hence, one would be 
able to cross check acceptances of the proton and neutron 
detectors. 

Note also that taking proton data from an independent 
run will potentially lead to another set of issues such as 
relative luminosity, the use of different beam energies, 
etc., which is likely to be on the level of 1%. 

In our analysis we neglect possible non-nucleonic com- 
ponents of the deuteron wave function such as kneaded 
(six-quark) components, etc. Current estimates put an 
upper limit on the probability of such components on 
the level of less than one percent, and they are expected 
to modify predominantly the high-momentum compo- 
nent of the deuteron wave function. Hence, the non- 
nucleonic components should give a very small (less than 
one percent) correction for the spectator momenta < 200 
MeV/c. Moreover, the experiments at the EIC looking 
for production of baryons such as A and N* in the spec- 
tator kinematics, x,l(A, N*) > 1/2, would allow to put 
a more stringent upper limit on (discover) non-nucleonic 
components of the deuteron wave function. 

In conclusion, we have demonstrated that a combined 
analysis of inclusive and semi-inclusive scattering off the 
deuteron coupled with a high resolution proton spectrom- 
eter will allow for the measurement of F^ at small x 
with the theoretical uncertainty of better than 0.5%. It 
would be a challenge to reduce the experimental system- 
atic errors to a comparable level. The measurement of 
the shape of the spectator spectrum would allow to de- 
termine nuclear shadowing in the deuteron with precision 
by far exceeding that possible in the inclusive measure- 
ments. 
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